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ABSTRACT 

A  review  of  procedures  for  designing  acceptance 
sampling  plans  and  process  control  schemes  according 
to  economic  criteria  is  given.  The  use  of  nonsta- 
tlsclcal  criteria  to  design  these  procedures  has 
bean  the  subject  of  considerable  research,  and  a 
number  of  different  process  models  and  procedures 
have  been  developed.  A  comparison  of  the  major  ap¬ 
proaches  to  model  formulation  and  design,  and  a  dis¬ 
cussion  of  the  practical  implementation  of  these 
techniques  is  the  primary  focus  of  this  paper. 

INTRODUCTION 

Acceptance  sampling  and  process  control  are  the 
primary  statistical  techniques  used  In  a  quality 
aesurahet  program.  These  techniques  have  found  wide 
application  in  industry,  particularly  in  manufactur¬ 
ing,  In  the  following  areas:  incoming  material 
inspection,  surveillance  of  the  production  process, 
estimation  of  lot  or  process  characteristics,  process 
capability  analysis,  and  finished  product  quality 
auditing. 

Traditionally,  acceptance  sampling  and  process 
Control  schemes  are  designed  with  respect  to  statis¬ 
tical  criteria.  For  example,  one  may  choose  the  sam¬ 
ple  size  and  acceptance  number  for  e  single- sampling 
plan  for  attributes  so  that  tha  operating  character¬ 
istic  curve  peases  through  (or  near)  two  points 
specifically  selected  to  give  certain  probabilities 
of  lot  rejection  at  specified  levels  of  lot  or  pro¬ 
cess  quality.  Similarly,  one  may  design  a  control 
chart  so  that  the  power  of  the  chart  to  detect  a 
particular  shift  in  process  quality  and  the  probabi¬ 
lity  of  false  alarms  are  equal  to  specified  values. 
While  tha  traditional  approach  often  produces  accep¬ 
table  results,  it  is  also  possible  to  design  accep¬ 
tance  sampling  and  process  control  schemas  with  re¬ 
spect  to  economic  criteria.  This  has  considerable 
intuit iva  appeal  since  there  has  been  an  increasing 
emphasis  on  quality  costs  in  racent  years.  Further¬ 
more,  the  aee  of  these  techniques  have  direct  econo¬ 
mic  consequences  in  that  one  is  balancing  tha  costs 
associated  with  sampling,  testing  and  process  sur¬ 
veillance  against  internal  and  external  failure 
costs.  Since  the  design  of  the  procedure  affects 
thaaa  coats,  it  is  logical  to  consider  this  design 
from  an  economic  viewpoint. 


During  the  last  20  years,  considerable  re¬ 
search  has  been  devoted  to  the  use  of  economic  cri¬ 
teria  to  design  acceptance  sampling  plana  and  pro¬ 
cess  control  schemes.  This  paper  gives  an  overview 
of  the  major  developments  and  approaches.  Some 
comparison  of  the  different  approaches  to  model 
formulation  and  a  discussion  of  the  practical  imple¬ 
mentation  of  these  techniques  is  also  given. 

ECONOMIC  MODELING  AND  PROCESS  CONTROL 

The  Shewhart  control  chart  la  probably  the  most 
widely  used  process  surveillance  device.  To  design 
a  control  chart,  we  must  choose  the  sample  size  (n) , 
the  control  limit  (k,  the  multiple  of  o) ,  and  the 
interval  between  aaaplea  (h  hours) •  Despite  the 
non-optimality  of  fixed  sample  sirs,  fixed  sampling 
Interval  procedures,  the  Shewhart  control  chart  has 
gained  widespread  use  because  of  its  flexibility, 
simplicity  of  administration  and  tha  additional  in¬ 
formation  about  process  performance  often  contained 
in  the  pattern  of  points  plotted  on  the  chart. 

Considerable  effort  has  been  devoted  to  deve¬ 
loping  economic  models  of  Shewhart  control  charts. 
These  models  usually  assume  that  tha  proceea  ia 
character iced  by  one  in-control  state  that  represents 
the  mean  of  the  quality  characteristic  whan  no 
assignable  causes  are  present,  and  s  *  1  out-of- 
control  states.  The  probability  model  chat  governs 
the  transit lona  between  these  s  +  1  stataa  is  called 
the  process  fsliura  mechanism. 

•  1 

Costs  and  Measures  of  Ef fectlvensaa 

Three  cetegories  of  costs  are  usually  consid¬ 
ered  in  the  development  of  economic  vodele  for 
Shewhart  control  charts;  tha  cost  of  sampling  and 
tasting,  usually  of  the  form  a^  +  S2n,  the  costa  of 
investigating  and  possibly  correcting  action  signals, 
i and  the  coat  of  producing  defective  items  (Internal 
and  external  failure  coats) .  These  coats  art  than 
combined  to  form  a  total  coat  par  unit  tima  function. 
Tha  general  approach  la  to  define  a  cycle  ea  tha 
length  of  tima  T  during  which  tha  proceaa  btglaa 
Optra ting  in  tha  in-control  state  and  eventually  re¬ 
turns  to  tha  ln-control  state  following  a  proceaa 
adjustment. 
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If  C  la  the  coat  incurred  during  a  cycle,  then 
the  expected  cost  per  unit  time  is 

EU)  -  E(C)/E(T)  (i) 

Since  E(C)  and  E(T)  are  functions  of  the  design  par* 
enters  of  the  control  chart,  numerical  optimization 
(direct  search)  techniques  car  be  employed  to  find 
the  parameter  values  that  minimize  the  expected  cost 
per  unit  time. 

Basic  Process  Models 

The  fundamental  process  model  often  used  in 
economic  design  of  the  Shewhart  control  c.*art  is  due 
to  (hmcan  (7j.  His  work  was  the  first  to  deal  with 
a  fully-economic  approach  to  control  chart  design 
and  to  explicitly  consider  the  optimization  problem. 
The  IXincan  model  assumes  a  single  assignable  cause 
(thus  a-1)  and  Chat  transitions  between  the  in- 
control  and  out-of-control  states  occur  according 
to  a  Poisson  process  with  intensity  X  occurrences 
per  unit  time.  Thus  the  length  of  time  the  process 
remains  in  control,  given  that  it  starts  in  that 
state.  Is  an  exponential  random  variable  with  mean 
1/A.  This  is  equivalent  to  a  process  failure  mech¬ 
anism  with  a  uniform  hazard  function. 

The  cycle  consists  of  four  periods:  (1)  the 
in-control  period,  (2)  the  out -of -control  period,  (3) 
the  time  to  taka  a  sample  and  interpret  the  results 
and  (4)  the  time  to  find  the  assignable  cause.  As¬ 
suming  that  the  process  continues  to  run  during 
searches,  the  expected  cycle  length  is 

K(T)  -  1/A  +  h/(X-0)  -  t  +  gn  +  D  (2) 

where  h  la  the  sampling  interval  (in  hours),  1  -  0 
Is  the  power  of  the  chart  for  detecting  a  specified 
out-of-control  stats,  n  la  the  sample  size,  g  is 
the  time  required  to  take  and  interpret  e  sample,  D 
Is  the  expected  time  to  find  an  assignable  cause  and 

1  (J+l)h  (j+l)h 

T  m  f  A(t-jh)dt  /  f  «“Xt  Adt 

Jb  JH 

!  -  [1  -  (l+Ah)*‘Xh]  /  A(l-._Xh)  (3) 

is  the  expected  time  of  occurrence  of  the  shift, 
given  that  it  occurs  betweeen  the  and  (j+l)4L 

samples. 

Let  Vo  end  Vj,  be  the  net  Income  per  hour  of 
operation  in  the  in-control  and  out-of-control 
states,  respectively.  The  costs  of  investigating 
real  and  falsa  alarms  are  a 3  and  el,  respectively. 
Siacs. there  era  E(I)/h  samples  per  cycle  and  oe'^V 
false  alarms  par  cycle,  on  the  average,  the 
expected  net  income  per  cycle  Is 

»<C)  -  V0(l/A)  +  (Sj+Sjft)  E(T)/h 
♦  Vx[h/(l-e)  -  x  +  gn  ♦  0] 

-  »j  -  «3««"Xh  /  a-«*Xh)  <*) 


Dividing  (4)  by  (2)  gives  the  expected  net  in¬ 
cone  per  hui.r,  say 

BCD  -  V0  -  E(L)  (5) 

where 

E(L)  -  Uj+a^/h 

+  Ja^[h/(l-0)  —  t  +  gn  ►  D] 

♦  *3  +  /  (1-e  )\ 

*  ll/A  +  h/(l-0)  -  t  +  gn  ♦  D]  (6) 

and  -  Vq  -  Is  the  penalty  coat  (par  hour)  re¬ 
sulting  from  production  in  the  out-of  control  state. 
Maximizing  E(I)  In  (5)  Is  equivalent  to  minimizing 
E(L)  in  (6).  There  are  three  decision  parameters, 
n,  k  and  h  (the  function  E(L)  depends  on  the  width 
of  the  control  limits  k  through  a  and  0);  and  nine 
user-supplied  constants  that  describe  the  production 
process,  a^,  S2»  *3,  *3*  *4* A,  g,  D,  and  $  the  mag¬ 
nitude  of  the  process  shift  (which  Is  Involved  In 
determining  o  and  0). 

While  this  basic  process  modal  could  be  uaad 
in  many  situations,  there  are  a  number  of  var let lone 
that  are  potentially  important.  The  first  of  theme 
Is  the  Incorporation  of  more  than  one  out -of -control 
state.  Duncan  [8]  has  generalised  the  cost  model 
above  to  a  situation  where  there  are  s  out -of -control 
atatee  (assignable  causes).  Knapp enberger  end  Cran- 
dage  [13]  also  present  a  model  capable  of  treating 
several  out-of -control  etataa.  The  Duncan  model 
appears  somewhat  more  realistic  in  that  it  employe 
different  search  costs  for  different  assignable 
causes,  while  the  Kneppenberger  and  Grandage  modal 
uses  an  average  or  composite  cost.  On  the  other 
hand,  thq  Kneppenberger  and  Grandage  modal  allow* 
continued  proc asa  deterioration  beyond  the  Initial 
shift,  while  the  Duncan  model  allows  only  one  or  two 
ehifte  to  occur.  Both  models  are  considerably  more 
sophisticated  In  terms  of  thslr  data  requirement a 
than  (6)  above.  Furthermore,  there  is  considerable 
evidence  Chat  a  tingle  assignable  cause  model  that 
matches  ths  true  multi-state  system  In  Important 
ways  will  be  a  satisfactory  approximation.  Specifi¬ 
cally,  if  we  let  the  single  out -of -control  state  be 
a  weighted  average  of  the  s  out-of-control  states 
with  the  weights  chosen  proportionally  to  the  proba¬ 
bility  of  occurrence  of  tha  etate,  excellent  results 
era  obtained.  In  practice,  there  is  probably  little 
reason  to  use  a  multiple  out-of -control  etate  model 
in  most  situations. 

Ths  assumption  of  s  constant  hasard  function 
for  the  process  failure  mechanism  (an  exponential 
distribution  for  the  elms  between  shifts)  is  criti¬ 
cal.  If  tba  occurrence  of  assignable  causes  can  ba 
thought  of  as  random  "ehocke"  disturbing  the  systma, 
than  tha  probability  of  a  ehlft  occurring  within  any 
email  interval  of  time  is  proportional  to  the  length 
of  the  Interval,  and  the  exponential  modal  is  prob¬ 
ably  appropriate.  However,  If  aaaigndble  esuaes 
occur  as  a  result  of  ths  cumulative  effects  of  vibra¬ 
tion  or  hast,  improper  set-up,  or  excessive  etressae 
during  start-up,  as  Is  ths  case  with  many  production 
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processes,  then  the  cxponrnCl.il  model  may  be  entirely 
Inappropriate.  Distr Ibutions  with  either  increasing 
or  decreasing  hazard  functions  may  be  more  reasonable. 
His spec If lcat Ion  of  this  aspect  of  the  process  model 
may  have  very  serious  economic  Impact.  Unfortunately, 
it  Is  not  straightforward  to  incorporate  other  process 
failure  mechanisms  in  (6).  A  discrete-t  line,  single¬ 
cause  model  that  allows  free  choice  of  the  process 
failure  mechanism  has  been  developed  by  Baker  [2], 
but  as  it  assumes  that  a  sample  is  taken  at  the  end 
of  each  period*  the  optimum  inter sample  interval  can¬ 
not  be  explicitly  determined. 

The  assumptions  in  (6)  that  the  process  contin¬ 
ues  in  operation  during  the  search  for  ar  assignable 
cause  and  that  the  co«t  of  repair  is  not  Included  are 
unrealistic  for  many  processes.  Fortunately,  (6)  may 
be  easily  modified  to  consider  these  possibilities. 

If  Dq  and  D2  are  the  expected  search  times  for  false 
and  real  alarms,  respectively,  and  if  &  is  the  cost 
of  repair,  then  asmimlng  that  the  process  Is  shut 
down  during  the  search  leads  to  an  expected  net  In¬ 
come  per  unit  time  of 

B(I)  -  -  E(L)  (7) 

where 

I(L)  -  j(«1+«2n)[l/X  +  h/ (1-6)  -  r]/h 

+  «3  ♦  A  +  [o«'Ah/(l-.'Ah)l(V0D0  +  *^) 

+  VpDjj  +  ^Ih/U-fi)  -  tj| 

♦  [i/x  ♦  h/u-M  -t  +  «d0*“x*7 

<l-*‘xh)  +  Dx]  (I) 

Tor  a  detailed  development  of  this  model*  see  H»t- 
gomery  [14).  Rote  that  given  a1#  a2,  *3,  a],  A,  a** 

1,  Vg,  Da*  D1#  and  the  magnitude  of  the  process  shift, 
($)  may  be  easily  optimised  for  n,  h,  end  k  using 
direct  search  methods. 

Wafer  Applications 

Suppose  that  the  quality  characteristic  of  in¬ 
terest  le  represented  by  the  perimeter  9,  end  that  x 
represents  the  sample  statistic  corresponding  to  9. 
Thera  ere  two  possible  values  for  9;  9  •  9 $  corres¬ 
ponding  to  the  in-control  state  sad  9  •  0^  which 
represents  the  out -of -control  state.  Then  the  pro¬ 
bability  of  e  false  elans  Is 

a  •  1  -  P(jLCL  gxi  DCL  I  •  *  9q)  C9) 

end  the  power  of  the  control  chert  Is 

1  -  $  -  1  -  TOCL  4  s  4  DCL  |  0  •  fj)  UO) 

where  ICt  sad  VCl  sad  tbg  lower  mod  upper  control 
limits  am  the  fhewhert  control  chert.  This  by  spe¬ 
cifying  the  quality  cherscc eristic,  the  sample  sta¬ 
tistic,  smd  tbs  relsvsmt  sampling  distribution*  the 
cost  models  of  tho  previous  section  could  he  applied 


to  any  type  of  Shewhart  control  chart.  For  example, 
if  the  quality  characteristic  of  interest  la  a  mea¬ 
surement  described  by  the  process  mean,  then  the 
x-chart  would  be  used;  thus  equations  (9)  end  (10) 
become 

a  *  24(-k) 
and 

1  -  8  «  +  *(-«yi-k) 

respectively*  where  9(r)  denotes  the  standard  normal 
distribution  function. 

The  major  applications  of  cost  models  such  as 
those  in  the  previous  section  have  been  to  the  x- 
chart,  the  x/R  chert  combination,  and  the  fraction 
defective  or  p-chart.  Some  work  has  also  focused  on 
procedures  for  the  simultaneous  control  of  several 
related  quality  characteristics,  and  the  treatment 
of  non-Shevrhart  control  charts  such  as  the  cumulative 
sum  control  chart  and  the  use  of  warning  limits  on 
the  x- chart .  Almost  all  reported  studies  assume  that 
the  process  failure  mechanism  is  exponential.  As 
noted  earlier*  this  is  a  critical  end  sometimes  un¬ 
warranted  assumption.  For  an  extensive  review  of  the 
literature  in  this  area*  see  Montgomery  [14]. 

Most  of  ths  papsrs  cited  In  fl4]  contain  numer¬ 
ical  studies  end  usually,  e  sensitivity  analysis. 

From  these  results.  It  is  possible  to  drew  certain 
general  conclusions  about  control  chert  design: 

1.  The  optimum  sample  else  is  largely  determined  by 
the  else  of  the  shift,  with  smaller  shifts  re¬ 
quiring  larger  swplas,  Smell  shifts  may  require 
very  large  samples,  possibly  as  large  aa  a*40  or 
more  on  an  x  chert. 

2.  Chang  as  in  e^  and  e2,  the  fixed  and  variable  costs 
of  sampling,  effect  ell  three  design  parameters. 
Increasing  tha  fixed  cost  a.  Increases  the  inter¬ 
val  between  samples  and  lamia  usually  to  slightly 
larger  samples*  while  increasing  the  variable  cost 
Sj  usually  results  in  Mil,  infrequent  samples* 
bat  narrow  control  limits. 

3.  The  penalty  coat  for  out -of -control  production 
•1  mainly  influences  the  interval  between  sam¬ 
ples  h.  Larger  values  of  a*  imply  smaller  opti¬ 
mum  values  of  h..  A  similar  sffact  la  observed 
by  increasing  X*  the  seen  number  of  process 
shifts  per  hour* 

4.  The  search  coats  aj  and  a«  mainly  affect  tha  con- 

!  trol  limits,  and  have  a  alight  affect  on  sample 

else.  Larger  values  of  03  sad  tj  result  in  wi- 
dar  control  limits  and  larger  temple  sleet,  aa 
increasing  star eh  testa  imply  that  fewer  false 
alarms  are  desirable. 

9.  The  optimum  economic  control  chart  design  is 
relatively  insensitive  to  estimates  of  the  coat 
parameters.  The  coat  surfaces  are  usually  flat 
naar  the  optimum,  end  are  usually  steeper  aaar 
tha  origin,  eo  that  it  la  bat  tar  to  overestimate 
the  control  chart  design  parameters  than  to  un¬ 
derestimate  them.  The  optlaum  design  la  moat 


sensitive  to  errors  In  ost treat the  in- 
control  and  ouc-of -control  states  (or  the  mag¬ 
nitude  of  th*  process  shift). 


Once  the  optical  n  and  k  are  determined,  the  Inter¬ 
val  between  samples  could  be  obtained  from 


This  last  finding  has  important  practical  impli¬ 
cations,  for  it  is  often  difficult  to  precisely  esti¬ 
mate  costs.  However,  the  mean  time  between  shifts 
and  the  magnitude  of  the  process  shifts  are  usually 
mora  easily  determined  from  process  performance 
data  or  from  the  engineer fs  knowledge  of  the  opera¬ 
ting  environment.  Note  that  one  should  exercise 
caution  in  using  arbitrarily  designed  control  charts 
and  empirical  "rules  of  thumb",  such  as  n» 5,  k*3  and 
h*l  for  the  x-chart.  In  some  cases,  particularly 
those  with  small  shifts  and  large  penalty  costs 
for  production  in  the  out-of-control  state,  signifi¬ 
cant  economic  penalties  may  result. 

A  Semi- Economic  Approach 

Instead  of  the  fully-economic  modeling  approach 
outlined  above,  it  may  be  desirable  in  some  cases  to 
use  a  semi-economic  scheme  that  blenda  both  economic 
and  statistical  criteria  directly.  This  could  be 
useful  in  situations  where  the  analyst  la  unwilling 
or  unable  to  estimate  all  of  the  fully-economic 
model**  cost  parameters,  or  where  he  wishes  a  com¬ 
promise  between  the  fully- economic  and  purely  statis¬ 
tical  approaches. 

One  such  semi-economic  scheme  would  be  to  find 
the  values  of  n,  k,  and  h  that  minimize  the  expected 
sampling  costs  per  cycle,  but  that  also  have  desir¬ 
able  average  run  lengths  for  the  control  chart  in 
both  the  In-control  and  out-of-control  states.  Con¬ 
sequently,  we  would  Ilka  to  choose  n,  k,  and  h  so  as 
to 

minimize  t  •  (a^+a^n)  E(T)/h  (11) 
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This  i9  an  approximation  for  the  optimum  sampling 
interval  h  given  n  and  k  suggested  by  Duncan  (7) 
and  Chiu  and  Wetherill  [A], 

These  formulations  of  the  optimal  control  chart 
design  problem  have  not  been  extensively  investi¬ 
gated.  It  would  be  of  interest  to  discover  hov  they 
compare  with  the  fully-economic  solutions,  and  with 
the  standard  designs  often  suggested  in  the  litera¬ 
ture. 


ECONOMIC  ASPECTS  OF  ACCEPTANCE  SAMPLING 


There  are  various  types  of  acceptance  sampling 
schemes;  single  sampling,  double  sampling,  sequen¬ 
tial  sampling,  continuous  sampling,  and  so  forth. 

The  quality  characteristic  being  Inspected  asy  be 
either  an  attribute  or  a  variable.  For  various 
reasons,  more  attention  has  been  given  to  attributes 
sampling,  and,  In  particular,  single  sampling  plans 
for  attributes  have  been  studied  extensively.  We 
now  give  some  of  the  key  results  pertaining  to  the 
economic  design  of  these  plans  and  briefly  survey 
some  of  the  other  related  work. 

An  attribute  single  sampling  plan  is  indexed  by 
three  numbers:  the  lot  size  K,  the  sample  size  n, 
and  the  acceptance  number  c.  If,  in  a  random  sam¬ 
ple  of  size  n,  c  or  fever  defectives  are  found,  the 
lot  is  accepted,  while  if  more  than  c  defectives  are 
found  the  lot  is  rejected.  The  probability  of  lot 
acceptance  is 

c  /N0\/N-Ne\  //N\ 

PC9)  •  l  (  4  ) (  n-d  )f  (n  )  (16) 


•here  ARX4  AHL^  are  the  average  run  lengths  in 
the  in-control  and  out -of -control  states,  respec¬ 
tively,  end  E(T)'  Is  given  by  equation  (2} .  Recall 
that  a  end  1-0  ere  functions  of  the  sample  size  n, 
the  control  limit  factor- k,  and  the  magnitude  of 
the  shift.  Rote  that  this  is  a  mult  ler  Iter  Ion  opti¬ 
misation  problem,  and  would  generally  he  harder  to 
solve  than  the  optimization  problems  usually  asso¬ 
ciated  with  Che  economic  design  of  control  charts. 
However,  it  avoids  explicit  estimation  of  the  search 
costs  S3  and  e^,  and  tha  penalty  cost  for  production 
la  the  out-of -control  state  04,  which  many  analyst t 
find  difficult. 

Another  possible  formulation  la  to  minimize  tha 
sampling  coat  on  a  per  sample  basis  while  simulta¬ 
neously  maximising  ARLg  end  minimising  ARLt .  Thus, 
<11)  becomes 

minimise  t"  •  a^  ♦  e^n  (1A) 

for  a  M,  end  <12}  end  Q3)  ere  unchanged.  This 
formulation  of  the  problem  Involves  only  n  and  k. 


where  6  Is  the  lot  fraction  defective.  The  plot  of 
?(6)  versus  6  for  0  s  ft  s  1  is  celled  the  operating 
characteristic  curve  (or  OC  curve)  of  the  plan. 
Traditionally,  we  design  acceptance  sampling  plans 
(l.e. ,  choose  n  and  c)  ao  that  P(6)  passes  through 
or  near  two  points,  such  as  the  familiar  producer's 
and  consumer's  risk  points  (AQL,a)  and  (LTPD,B) . 

Whan  N  Is  large  relative  to  n,  the  hyper geometric 
distribution  in  (16)  can  be  replaced  by  the  binomial. 

Prior  Plat rlbut lone  for  Attribute  Sampling 

Lot  fraction  defective  8  la  a  function  of  two 
sources  of  variability;  the  variability  of  p  tha 
process  fraction  defective,  end  the  variability  of 
6  about  p.  It  la  usually  convenient  to  assume  that 
lot  quality  has  e  mixed  binomial  distribution;  chat 
la,  each  lot  Is  produced  by  a  production  process 
that  is  in-control  at  the  level  p,  but  p  varies  from 
lot  to  lot  according  to  a  probability  distribution 
f(p).  The  distribution  f(p)  la  often  called  the 
prior  distribution  for  p  or  the  process  curve.  It 
la  extreeely  important  to  note  that  If  the  process 
quality  la  stable  such  that  f(p)  •  1  whan  p  •  p(j 


and  f(p)  •  0  when  p  ^  pq,  th«n  '•here  is  no  need  for 
sampling  (this  result  is  called  Mood's  theorem;  see 
*>od  [15]). 


where  a  is  d  constant  proportional  to  the  variable 
cost  of  sampling  and 


While  there  are  many  possible  choices  tor  the 
prior  distribution  for  p,  some  of  the  more  important 
are  the  following:  the  continuous  beta  distribution. 


?<p>  *  I  (d)pda-p)n"d 
d-0  v  ' 

is  the  probability  of  accepting  a  lot  of  quality  p.  ’  « 
Minimizing  L  with  respect  to  n  and  c  will  produce 
the  optimal  sampling  plan.  This  is  sometimes  celled 
the  Bayesian  approach  to  designing  a  sampling  plan. 


the  discrete  two-point  binomial, 

PCp-pj)  *  fi»  1*1,2,  (18) 

and  the  normal-generated  distribution 

f(p)  -  cT1  expl^p2  -  *s(u-» )2/o2]  (19) 

where  p  *  ♦(-u) .  In  practice,  it  Is  important  to 
know  how  accurately  the  prior  distribution  must  be 
epee if led.  Generally,  the  analyst  does  not  possess 
sufficient  information  about  the  process  to  specify 
the  nr lor  with  great  confidence.  Fortunately,  most 
results  indicate  that  precise  specification  of  the 
prior  la  not  critical,  provided  that  a  reasonable 
distribution  la  chosen.  Continuous  prior  distribu¬ 
tion*  are  generally  thought  to  be  more  appropriate 
than  discrete  ones,  end  the  beta  distribution  (17) 
has  been  used  extensively.  However,  when  the  under¬ 
lying  quality  characteristic  Is  a  continuous  variable 
that  la  normally  distributed  within  each  lot  and  the 
mean*  of  this  quality  charactt fistic  also  has  a  nor¬ 
mal  prior,  then  tha  prior  distribution  for  p  has 
tha  form  (19) .  The  beta  and  normal -genera tad  dis¬ 
tributions  can  have  very  different  shapes  for  the 
same  mean  and  variance,  and  so  slgaif leant  differ¬ 
ences  In  the  optimal  sampling  plans  may  result.  For 
further  discussion  of  prior  distributions,  see 
Vetherlli  and  Chiu  [21]. 


The  Major  A 
Attributes 


prose has  to  Singls  Sampling 


Given  a  suitable  prior  distribution  and  e  set 
of  costs  or  losses  associated  with  sampling  plan 
operation,  it  la  desirable  to  choose  the  sampling 
plan  parameters  that  minimise  the  total  cost.  Per¬ 
haps  tha  most  widely  used  and  detailed  model  is 
that  of  Guthrie  and  Johns  [9],  A  simplified  form 
of  this  model  is  also  presented  by  Raid  [XI].  The 
model  la  a  linear  cost  model.  All  linear  cost  for¬ 
mulations  lead  to  tha  same  expected  loss  function. 

Ha Id  utilises  the  concept  of  break-even  quality  p^, 
a  fraction  defective  value  at  which  It  le  Just  as 
costly  to  sccept  as  to  reject  the  lot.  The  expected 
lose  per  lot  Is 


‘■-Hi  lpt*J>)  f  1— »Cpl  ]t  lp)dp 
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The  major  focus  of  the  research  in  this  area  has 
started  with  the  Guthrie  and  Johns  model.  Held  has 
been  a  major  contributor  in  the  field,  along  with 
some  of  his  coworkers.  The  major  emphasis  has  been 
on  finding  asymptotic  relationships  between  n  and  c 
for  various  process  curves,  and  in  producing  tables 
suitable  for  use  by  professional  practitioners.  Hald 
[10]  gives  a  number  of  significant  findings  in  his 
1960  paper.  One  part  of  the  paper  Investigates  the 
compound  hypergeometric  distribution;  that  is,  the 
probability  distribution  of  the  number  of  defectives 
d  in  a  random  sample  of  sice  n  given  a  prior  distri¬ 
bution.  The  second  part  of  the  paper  is  directed 
towards  actually  finding  optimum  sampling  plana  for 
rectangular,  beta,  and  double  binomial  priors.  A 
general  solution  la  given  assuming  the  linear  cost 
model  and  inequalities  are  given  for  n  and  c.  In 
1965,  Hald  [111  provided  tables  for  the  double  bino¬ 
mial  prior  and  In  1966,  Hald  [12]  provided  tables 
fox  the  beta  prior,  along  with  asymptotic  relation¬ 
ships  between  n  and  c,  and  between  N  and  n. 

Despite  the  significance  of  Held’s  work,  the 
tables  he  has  provided  are  often  difficult  to  use 
because  of  the  large  amount  of  information  required. 
An  alternative  approach  consists  of  formulating  an 
appropriate  cost  model  and  optimising  it  for  a  spe¬ 
cific  problem  using  direct  search  methods.  Consi¬ 
derable  work  in  thla  area  has  been  done  by  G.K.  Ben¬ 
nett,  K.B.  Case,  and  J.V.  Schmidt  and  their  students. 
For  example.  Chair  1972  and  1975  papers  [3] [17]  deve¬ 
lop  economic  models  for  single  sampling  plans  for 
dealing  simultaneously  with  multiple  attributes.  The 
coat  models  consist  of  a  component  representing  In¬ 
spection  coats,  a  component  representing  the  expec¬ 
ted  coat  of  lot  rejection,  and  a  component  represen¬ 
ting  the  Expected  coat  of  lot  acceptance.  Fat  tern 
search  is  used  for  model  optimisation.  In  general, 
direct  search  methods  are  a  very  effective  approach 
for  determining  economically  optimal  acceptance 
sampling  plans. 

The  economic  impact  of  the  disposition  policy  for 
rejected  lots  has  also  been  investigated.  Vhlle 
there  are  a  numbar  of  posslbla  lot  disposition  pol- 
iclaa,  the  tiro  cases  that  have  bean  Investigated  moat 
extensively  ere  where* rejected  lots  are  either 
scrapped  or  screened  (100  percent  Inspection).  For 
work  on  lot  disposition  policies,  see  [18]  end  [20]. 

Other  Work 

There  have  been  many  other  studies  devoted  to 
the  economic  design  of  acceptance  sampling  plans; 
for  Instance,  a  1975  survey  paper  by  Wat her 11 l  and 
Chiu  [21]  cites  25)  references.  Vhlle  most  of  the 
work  focuses  on  singls  sampling  for  attributes,  there 


has  been  ecxae  research  on  the  economic  design  of 
•ore  sophisticated  sampling  plans*  Foe  example, 
Stewart,  Montgomery  and  He ikes  [20]  describe  a  proce¬ 
dure  for  the  selection  of  double  sampling  plans  for 
attributes  based  on  prior  distributions  and  costs. 
Models  are  presented  for  the  cases  where  rejected 
lots  are  either  screened  or  scrapped*  They  note  chat 
there  is  often  little  difference  in  cost  between 
economically  optimal  single  and  double-sampling  plans. 
However,  when  sampling  cost 3  are  large,  double  sam¬ 
pling  plans  have  much  to  offer.  They  also  observe 
that  arbitrary  double- sampling  plans,  such  as  those 
in  MIL  STD  105D,  may  be  very  far  from  economically 
optimum. 

The  effect  of  inspection  error  on  sampling 
plan  design  has  also  received  considerable  attention. 
If  an  inspector  misclasslf let  good  and  bad  items  with 
constant  probabilities,  the  effect  is  to  translate 
the  0C  curve  of  the  sampling  plan,  so  that  the  actual 
or  effective  0C  curve  is  somewhat  different  from  the 
nominal  or  advertised  0C  curve.  If  the  probabilities 
of  mlsclasslf icatlon  are  known,  then  one  may  directly 
incorporate  thla  information  into  the  economic  design 
of  the  sampling  plan.  Generally,  the  presence  of 
inspection  errors  implies  that  larger  samples  are 
necessary.  A  good  review  of  the  literature  in  this 
area  Is  in  Dorris  and  Foote  [6], 

Very  little  ettention  has  been  given  to  accep¬ 
tance  sampling  by  variables.  Variables  sampling  is 
not  as  widely  used  In  practice  as  attributes  sampling 
for  several  reasons: 

1*  Variables  measurement  is  often  more  expensive 
(difficult)  than  attributes  measurement. 

2.  A  separate  plan  must  be  used  for  asch  quality 
characteristic. 

3.  The  estimation  of  fraction  defective  assumes  a 
normally  distributed  quality  characteristic.  If 

.  this  assumption  is  violated,  the  tail  areaa  may 

!  be  dramatically  affected  and  the  resulting  frac¬ 

tion  defective  estimate  grossly  in  error. 

4.  It  la  possible  to  reject  a  lot  without  actually 
finding  any  defectives,  and  this  often  upsets 
both  producers  and  consumers. 

However,  remember  that  variables  sampling  can  greatly 
reduce  the  required  sample  site,  and  that  It  docs 
generally  provide  betteT  information  about  the  lot 
or  process  quality.  For  work  on  the  economic  design 
of  variables  sampling  plan,  see  Allor,  Schmidt  and 
Bennett  [1],  Schmidt,  Bennett  and  Case  ilB],  and 
Scla&ldt,  Case  and  Bennett  [19].  Reference  [1]  deals 
with  the  situation  where  the  quality  characteristics 
are  a  mixture  of  attributes  and  variables.  In  all 
of  these  studies,  ths  approach  taken  is  to  formulate 
a  cost  model  end  optimise  it  via  direct  search 
methods. 


The  last  20  years  have  seen  the  development  of 
aumetoue  techniques  for  the  design  of  process  control 
end  acceptance  samp ling  schemes  based  on  economic 
considerations.  However,  the  indication  is  that  very 


very  few  practitioners  have  implemented  any  of  these 
techniques  (for  example,  see  the  surveys  by  Sanlgs 
and  Shirland  [16]  and  Chiu  and  Wetherlll  [5]).  This 
is  surprising,  as  most  quality  assurance  managers 
claim  that  cost  reduction  and  increased  productivity 
Is  a  major  objective  of 'their  function.  In  many 
cases,  an  experienced  engineer  could  design  an  appro¬ 
priate  technique,  perhaps  even  one  that  is  nearly  • 
economically  optimal,  but  the  use  of  a  formal  econo¬ 
mic  model  to  assist  the  analyst  Is  a  much  more  pre¬ 
cise  approach,  leaving  leas  to  judgment.  Further¬ 
more,  there  is  often  a  significant  economic  penalty 
associated  with  the  "standard",  judgment  designs  so 
frequently  used  in  practice. 

The  implementation  of  these  techniques  requires 
a  computer  program  of  the  cost  model  and  the  optimi¬ 
zation  procedure.  The  lack  of  availability  of  suit¬ 
able  computer  software  has  certainly  slowed  practical 
implementation.  This  is  an  Important  gap  between 
theory  and  practice  that  must  be  filled  before  the 
economic  design  of  quality  assurance  techniques  will 
become  widespread. 

Many  practitioners  are  reluctant  to  use  these 
techniques  because  of  the  difficulty  In  estimating 
costs,  prior  distributions,  and  other  model  parame¬ 
ters.  Fortunately,  coats  and  prior  distributions  do 
not  have  to  be  estimated  with  high  precision,  although 
some  other  parameters,  such  as  the  magnitude  of  the 
process  shift,  require  more  careful  determination. 
Sensitivity  analysis  of  the  specific  model  could 
help  the  analyst  discover  which  parameters  are  cri¬ 
tical  in  bis  specific  application.  The  availability 
of  efficient,  interactlva  computer  software  would  be 
of  significant  value  in  this  respect. 

.  He  will  illustrate  the  use  of  a  simple,  inter¬ 
active  computer  program  for  the  optimal  economic 
design  of  an  x-chart .  The  program  assumes  that  tha 
Duncan  single-cause  process  model  la  appropriate, 
and  requests  the  user  to  Input  values  of  a^,  *2  *  *3* 
al,  *4,  1|  6*  g,  and  D,  The  optimal  control  limit 
width  k,  interval  between  samples,  and  minimum  cost 
per  unit  time,  are  calculated  and  displayed  for  a 
range  of  sample  sizes,  along  with  the  o-risk  end 
power  of  tne  control  chart  for  each  (n,k,h)  combina¬ 
tion.  The  economically  optimal  control  chart  design 
may  be  found  by  inspection  of  the  cost  function  val¬ 
ues  to  find  the  minimum.  The  output  provided  enables 
the  analyst  to  determine  the  sensitivity  of  the  cost 
surface  In  the  vicinity  of  the  optimum.  The  program 
Is  written  In  FORTRAN  for  a  CDC  CYBER- 7 4  computer, 
and  involves  leaa  than  100  lines  of  code.  Execution 
times  are  typically  one  second  or  less.  Further 
details  of  the  program  are  available  from  the  author. 

Consider  a  manufacturer  of  non-returnable 
glass  bottles  for  packaging  a  carbonated  soft  drink 
beverage.  The  wall  thicknase  of  the  bottles  lr  an 
important  quality  char actar  1st lc.  If  the  wall  1* 
too  thin,  internal  pres aura  generated  during  filling 
will  causa  ths  bottle  to  burst.  The  manufacturer 
has  been  using  x  and  R  charts  with  n*3,  k*3 ,  and 
taking  samples  every  h»2  hours  to  control  the  pto- 
caaa.  He  wishes  to  compare  this  with  an  aconcmlcally- 
opt Intel  design  and  satinet*  the  savings. 

Baa  ad  on  an  analysis  of  ths  quality  control 


technicians*  salaries,  the  costs  of  the  test  equtp- 
neot,  and  the  length  of  time  requl^d,  it  i« 
ted  that  the  fixed  cost  of  taking  a  sample  id  61.00, 
and  the  variable  cost  of  sampling  Is.  approximately 
$0.10  per  bottle.  It  takes  about  one  minute  (0.0167 
hours)  to  measure  and  record  the  wall  thickness  of  a 
bottle. 

The  process  is  subject  to  several  types  of 
assignable  causes  resulting  in  large  shifts,  typi¬ 
cally  of  about  two  standard  deviations.  (This  is 
why  the  "standard"  x-chart  design  n-5,  k*3  has  been 
used  -  it  is  known  to  be  reasonably  effective  for 
large  shifts).  The  previous  six  months  of  line  per¬ 
formance  data  Indicate  that  the  mean  time  between 
process  breakdowns  is  about  20  hours.  Thus  \»0.05  is 
the  parameter  of  the  exponential  distribution.  The 
average  time  to  investigate  an  action  signal  is  one 
hour.  Real  alarms  Incur  a  $25.00  search  cost,  on 
the  average,  vhlla  falsa  alarms,  which  are  more  dif¬ 
ficult  to  check  out,  incur  a  search  cost  of  $50.00. 

The  soft  drink  bottler  to  whom  the  bottles  are 
sold  has  a  policy  of  backcharging  the  bottle  manu¬ 
facturer  for  the  costs  of  cleanup  and  lost  production 
when  an  excessive  number  of  defective  bottles  burst 
during  filling.  The  past  six  months'  experience 
indicates  that  $100  per  hour  la  a  reasonable  estimate 
of  the  penalty  cost  of  operating  in  the  out-of -control 
state. 


Figure  1  snows  the  actual  computer  input  and  out¬ 
put  for  this  example.  Note  that  the  optimum  design 
has  n«5,  k*2.$y,  h*0.76  hours  (about  45  minutes),  and 
•  “inimum  cost  of  $10.38  per  hn^r.  Thus  the  implica¬ 
tion  is  that  the  manuf'^urer’s  existing  control  chart 
w’sign  is  good  w**h  respect  to  n  and  k,  but  that  he 
need**  to  samp,  more  frequently.  The  actual  cost  of 
his  cur*,-**#-  procedure  Is  $12.05  per  hour,  or  about  ’a 
16  percent  penalty  cost.  While  the  savings  per  hour 
seem  small,  this  is  a  continuous  production  process 
operating  three  shifts  per  day,  and  the  annual  sav¬ 
ings  are  over  $14,000.00. 

»  . 

After  looking  at  the  optimal  x-chart  design,  the 
bottle  manufacturer  suspects  that  he  may  have  incor¬ 
rectly  estimated  the  penalty  cost  of  out-of-control  . 
production  (a^),  and,  at  worst,  he  may  have  underes¬ 
timated  this  parameter  by  50  percent.  Therefore,  he 
reruns  the  program  with  a4«$150.00  to  investigate 
the  effect  of  mlsspeclfylng  this  parameter.  The  com¬ 
puter  output ,  shown  in  Figure  2,  indicates  that  the 
optimum  design  is  now  n-5,  V-2.99,  h-0.62,  with  a 
minimum  cost  of  $13.88  per  hour.  Note  that  the  pri¬ 
mary  effect  of  increasing  a$  by  50  percent  is  to  de¬ 
crease  the  interval  between  samples  from  45  minutes 
to  37  minutes.  The  program  could  be  used  in  this 
manner  to  quickly  and  easily  investigate  the  effects 
of  errors  in  specifying  any  of  the  model  parameters. 
Based  on  this  analysis,  the  manufacturer  elects  to 
adopt  a  45  minute  sampling  interval,  because  it  is 
administratively  simple  and  it  will  not  require  any 
additional  quality  control  technicians. 
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Figure  1.  Sample  Computer  Oitput 
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Figure  2.  Sample  Computer  (Xitput 
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